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Abstract 

A new formulation of the spectral energy budget of kinetic and available potential energies of the atmosphere 
is derived, with spherical harmonics as base functions. Compared to previous formulations, there are three main 
improvements: (i) the topography is taken into account, (ii) the exact three-dimensional advection terms are considered 
and (iii) the vertical flux is separated from the energy transfer between different spherical harmonics. Using this 
formulation, results from two different high resolution GCMs are analyzed: the AFES T639L24 and the ECMWF 
IFS T1279L9I . The spectral fluxes show that the AFES, which reproduces realistic horizontal spectra with a fc _5//3 
inertial range at the mesoscales, simulates a strong downscale energy cascade. In contrast, neither the k~ 5 ^ 3 vertically 
integrated spectra nor the downscale energy cascade are produced by the ECMWF IFS. 



Introduction 



i 

C3 The atmospheric horizontal spectra of velocity components 
c/^and temperat ure show a robust k~ 5 ^ 3 r ange at mesoscales 
# (10-500 km) (jNastrom and GaeeL Il985h . This power law 
c/)and the corresponding one for the horizontal second or- 
£^>ier structure functions (Y 2 / 3 ) are obtained from signals 
*^ measured in the troposphe re and the stratosphere, ove r 



i__i)oth land and sea (see e.g., Frehlich and Sharman . 201dh 
It still remains a challenge to reproduce these results in 
^ simulations. Some genera l circu l ation models (GCMs) 



(e.g. iKoshvk and Hamilton! . l200ll ; lHamilton et all . l2008h 



and mesoscale numeri cal weather prediction (NWP) mod- 
^sO els (|Skamarockl 120041) reproduce quite realistic mesoscale 
0> spectra. Other GCMs, as for example ECMWF's weather 
1— ( prediction model Integrated Forecast System, produce 
mesoscale spectra significantly steeper and with smaller 
magnitude than the measure d ones, even w ith relatively 
. .high resolution versions (e.g., IShuttd . 120051) . The inabil- 
J>"rty of some GCMs to simulate realistic mesoscale spectra 
k>( mu st have impo rtant consequences in terms of predictabil- 
*^jrty ( Vallisl. 2006), dis persiveness of ensemble prediction sys- 
terns (jPalmerL I2OO 1 ) and, evidently, mesoscale NWP. 

Even though one can now simulate realistic mesoscale 
spectra, it is still unclear what physical mechanisms pro- 
duce them. Theoretically, the only convincing explana- 
tion of the fc -5 / 3 power law is the hypothesis that it 
is produced by an (upscale or downscale) energy cas- 
cade with a constant energy flux through the scales. 
Therefore, most of the different theories proposed are 
based on the hypothesis of an energy ca scade: up- 
scale cascad e due to 2D-stratified turbulence ( GaeeL 19791 : 



Lilly], 1198311 or downscale cascade due to internal gravity 



waves ( Dewan and Goodl . 119861: ISmith et all . Il987h . quasi- 



geostrophic dynamics ( Tung and Orlando . 20031), surf ace- 
quasigeostrophic dynamics ( Tulloch and Smith!. 20091 ) or 
3D-strongly stratified turbulence ( LindborgLl2006 ). 



The principle of energy conservation strongly constrains 
the dynamics of the atmosphere. In or der to explain the 



Lorend (|l955l) de- 



maintenance of the general circulation, 
veloped the concept of available potential energy (APE) 
and derived an approximate expression proportional to 
the variance of the temperature fluctuation. His work 
has laid th e foundations for seve r al studies invest i gating 



APE (e.g.. | Boerl Il989l: IShepherdl. fl99l : ISiegmundl . Tl994 : 
Molemaker and McWilliamsl . 12010 ) and the atmospheri 



nc 



energy budget through diagnostics of d ata from global me- 
teoro l ogical analysis and GCM s (e.g., Boer and Lambert . 
2008 1 ; ISteinheimer et all . 120081 ). Due to the multiscalc na- 



ture of atmospheric motions, spectral analysis c an revea l 
valua ble pieces of information from the data (jFiOrtoftl . 



1953) and have therefore become a standard method for 
diagnostics. However, drawbacks of different used formu- 
lations of the spectral energy budget severely limit the 
results. First of all, the attention has been focused on 
the budget of kinetic energy (KE) so that in many stud- 
ies the APE budget is not considered. Most studies in- 
vestigate only the budget integrated over the total height 
of the atmosphere and thus the vertical fluxes of energy 
are not computed. Another very important limitation is 
tha t most s t udies are based on the formulation proposed 
by iFiOrtoft ( 1953 ). in which only the p urely horizontal 



Boer and Shepherd. 1983; 


SheDherdl. 119871 iBoerl. 


1994; 


Straus and Ditlevsenl. 1999; 


Burgess et al.. 2012). This av- 



1976 



proximation is justified for the very large scales of the at- 
mosphere for which the divergence is indeed very small. 
However, the approximation may lead to large errors at 



* Corresponding author address: Pierre Augier, Linne Flow Centre, KTH Mechanics, Stockholm, Sweden. 
E-mail: pierre.augier@mech.kth.se 



Preprint submitted to J. Atmos. Sci. 



October 15, 2012 



1 



the mesoscales. Atmospheric measurements show that 
divergent and r otational spectra are of the same order 
( Lindborg] . 2007 ) and atmospheric simulations produce di- 
vergent spectra of the same order of magnitude as rota- 
tional spectra at the mesoscales (e.g. . Hamilton et al.l . 12008: 



Burgess et all 12012 ). Moreover, these results are consis- 



tent with theoretical results showing that strongly strati- 
fied and weakly rotating flows tend to evolve toward states 
in which the divergent and ro tational components are of 
the same order of magnitude (Billant and Choma3 . 2001 



Lindborg and BrethouwerL 120071 ). Therefore the spectral 
energy budget formulations based on the two-dimensional 
vorticity equation can not capture the dynamics at the 
mesoscales. 

There is no theoretical obstacle in considering the exact 
three-dimensional advection includi ng both ro t ation al and 
divergent components of the flow. Lambert! ( 1984 ) have 
developed a formulation of the spectral energy budget con- 
sidering both KE and APE and taking into account the 
exact advection term. However, the diagnostics was inte- 
grated over the total height of the atmosphere so vertica l 
fluxes were not considered. Koshvk and Hamilton! ( 2001 ) 
performed a diagnostic of the equation for the KE spectrum 
(the APE budget was not investigated). The exact advec- 
tion was computed but the vertical flux was not separated 
from the energy transfer between spherical harmonics, so 
it wa s impossible to define spectral fluxes in a conservative 
way. iKoshvk and Hamilton! (|2001l ) separated the spectral 
pressure term into adiabatic conversion and vertical flux. 
However, the separation was only approximate. Moreover, 
as in all previous studies on the spect ral energy budget, 
the to pography was neglected. Recently. iBrune and Becker 
( 2012! ) have investigated the effect o f the vertical resolution 
in a mechanistic GCM. Just as in Koshvk and Hamilton 



(|2001l ). all the terms in the kinetic spectral energy equa- 
tion were computed at different pressure levels and it was 
demonstrated that they balance each other. However, spec- 
tral fluxes were defined in a non-conservative way and ac- 
tually also included vertical fluxes. 

In order to investigate the energetics of the mesoscales 
simulated by GCMs, it would be desirable to formulate the 
spectral energy budget considering both KE and APE, tak- 
ing the topography into account and making an exact sep- 
aration of the advection terms into spectral transfer and 
vertical flux and a corresponding separation of the pres- 
sure term into adiabatic conversion, vertical flux and spec- 
tral transfer. A formulation meeting these requirements is 
derived in section [31 In section [31 results from two high 
resolution GCMs are analyzed. 

2 Formulation of the spectral energy bud- 
get 

2.1 Governing equations in p- coordinates 

The analysis is performed in pressure-coordinates in which 
variables are functions of time, longitude A, latitude ip (the 
horizontal coordinates are denoted by x/j) and pressure. 
The main advantage of the p-coordinates is that mass con- 
servation can be expressed in the same way as for an in- 



compressible fluid: V • v = 0, where V = (Vh,d p ) is 
the gradient operator and v = (u, ui) is the total veloc- 
ity, with the horizontal gradient operator, u the hor- 
izontal velocity and u> — D t p the pressure velocity (D t 
is the material derivative). The hydrostatic equation is 
d p & = —a = —RT/p, where $ is the geopotential and a 
the volume per unit mass. 

The potential temperature is given by 9 = A(p)T, with 
Hp) = (Po/p) x , Po = 1 bar, x = R/cp - 2 /7- For sim- 
plicity, the corrections related to the vapor content are 
neglected. However, the latent heat release is taken into 
account through the rate of production of internal energy 
by heating Q. The main drawback of the p-coordinates is 
the complication related to the lower boundary condition 
(the topography pierces the lower pressure le vels) . It can 
be overcome with the formalism developed bv iBoer (1982) 
(see appendix) in which the dynamical equations can be 
written as 

d t u = - v • Vu f{<p)e z A u — pV h § + 0D u (u), (1) 
d t & =-v-W- ud p (9) r + Q e - I3d t (0) r + PD g (0), (2) 

where f(<p) is the Coriolis parameter, Qg = A{p)Q / 'c p , 
u = /3(x/j,p)u, with (3(xih,p) equal to one above the surface 
and to zero below. The potential temperature fluctuation 
§' is defined as 9' = 0-/3(9) r , where (9) r = (/30)/(/3) is the 
representative mean, i.e. the mean over regions above the 
surface (the brackets (•) denote the mean over a pressure 
level). D u (u) and Dg(#) are diffusion terms. 

2. 2 Kinetic and available potential energy forms 
lLorenzl (fl955l) showed that the sum of the globally inte- 



grated kinetic energy (KE) and available potential energy 
(APE) is approximately conserved. The mean energies per 
unit mass are E K {p) = (|u| 2 )/2 and E A (p) = j(p)(9' 2 ) /2, 
with j(p) = R/(—A(p)pdp(9) r ). The energy budget can be 
written as 

d t E K {p) = C(p)+d p F Kt (p)-D K (p), (3) 

d t E A (p) = G(p)-C( P )+d p F At (p)-D A (p) 

+J(P), (4) 

where G(p) = r ){p){9'Q l g ) is the forcing by heating (differ- 
ential heating at very large scales and latent heat release), 
C(p) = —{Qa) is the conversion of APE to KE, F A f(j>) = 
1 {p)(tu9' 2 )/2 and F Kt {p) = -(^|u| 2 /2) - + (d tPs <Z> s p) 
are the vertical fluxes and Dk{p) and D A {p) are diffu- 
sion terms. The last term of Q, J(p) = {d p ^)(uj9' 2 ) /2 - 
(/3 2 )(w) r (a) r , corresponds to adi abatic pro c esses which do 
not conserve the Lorenz APE. ISiegmundl ( 1994 ) showed 
that the globally integrated exact APE and Lorenz APE 
differ by less than 3%. 



2.3 Spectral analysis based on spherical harmonic trans- 
form 

Each function defined on the sphere can be expanded as 
a sum of spherical harmonics functions Y; m (x^), the nor- 
malized eigenfunctions of the horizontal Laplacian operator 
on the sphere: \V h \ 2 Y im = l(l + l)Y lrn /a 2 and (Y^Yim) = 
o~ii'0~mm' , where the star denotes the complex conjugate and 
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/ a nd m are th e total and zonal wavenumbers (for details, 



see 



Boeii 119831 ) . The velocity is written as 



= ^ ^ ui m (p)Yi 



(5) 



The nonlinear KE and APE spectral transfer terms can 
be expressed as 

Tl m \p) = U{ - (uL ■ [u ■ V h u] lm + u; m ■ [uV, • u] Jwi /2) 



l>0 -l<m<l 



((d p u 



m]l m - ui m ■ [cvd p u]i m ) /2}, 

(11) 



and the other variables are writtah the corresponding T { j m] (p) = 5R{ - 7 (^[u • V^'] /m + 0£J0'V h • u], m /2 



way. The spectral KE function E^ n> (p) — |u ;m (p)| 2 /2 is 
defined so as the mean kinetic energy at a pressure surface 
is given by E K {p) — J2i m ^k^ip)- F° r APE, the equiva- 
lent requirement gives E^\p) = r y(p)\d' hn (p)\ 2 /2. 

The spectral energy budget is derived by multiplying 
equation (fTJ) by u* m and equation @ by "l{p)0[* m and then 
taking the real part of the resulting equations. In order 
to separate the pressure term into conversion, transfer and 
vertical flux, we use the continuity equation and the hydro- 
static equation to rewrite it as 

- fi L • [/3V ft $]; m = - tf m a lm - d p (w* lm $>i m ) 

-tf m [5<S> s ]i m - [d tPs 8]* lm $>i m 

+ u^ ro • [$ s V h p s 5] lm , (6) 

In order to separate the advection terms into transfer and 
vertical flux we rewrite them as 

- VL* lm ■ [u • V h u]i m - VL* lm ■ [uVft • u] /m /2 
+ ({d p U* m ) ■ Hi™ - U*m ' [ud p u]lm)/2 

- 9 p (u* m ■ [wu]im)/2, (7) 
- 7 (p)t[vV4 m = 

- 7 CJu ■ V h 9'] lm - !~eT m [0'V h ■ u] lm /2 

+ 7 ((^)M'W - e'c m [ud p e'] lm )/2 
-a p ( 7 c„M'] ;m )/2 

+ (d p ~,)e' l * m ["0'}lm/2, (8) 

where we have used the continuity equation. Neglecting 
the non-conservative term in ([5]) corresponding to the last 
term in Q , the spectral energy budget can now be written 
as 



\P) 



C [lm] ( P ) 



rp[ln 
1 K 



Kp) 



-d v F^\p) 



sl m \ P ) 



dT\p) 



d t E^\p) 



G^ lm \p)-C^ m \p) + Tl m \p) 

-J/rn 



L^(p) 
(9) 



-dX? ] (p) 



where T^ m \p) and T^ n \p) are spectral transfer terms 
due to nonlinear interactions, L\ lm \p) is a spectral trans- 
fer term arising from the Coriolis and pressure terms, and 
S K m ^ (p) is a surface term. Each of the other terms cor- 
responds to a term in (J5H3]). In contrast to the previous 
formulations, here, the APE budget is included, the topog- 
raphy is taken into account, the exact three-dimensional 
advection is considered and the vertical flux terms and the 
horizontal spectral transfer terms are exactly separated. 



(12) 

where 3? denotes the real part. The transfer term arising 
from the Coriolis and pressure terms is 

L^(p) -(u? ro -[V h $] Jni + [V fc -u]r m $ Iw ,) 

-u; m .[/^)e z Au] im }. (13) 

Each transfer term, as each of the sums of the terms of 
each line of (fTTjl . ([12"]) and (TIB"]) , is exactly conservative, 
meaning that it only redistributes energy among the dif- 
ferent spherical harmonics at a pressure level. For exam- 
ple, the sum over all spherical harmonics of the first and 
second lines of (fTTjl are equal to (—'Vf l {u\u\ 2 )) = and 
((d p u)uju — uu(d p u)) = 0, respectively. The terms in (fT5|) 
have been gathered together because in plane geometry, 
where a Fourier decomposition is used, the corresponding 
terms are zero. Note also that for purely geostrophic mo- 
tion, these terms cancel each other (except at the surface) 
and L^(p) ~ 0. 

The diabatic term, the conversion term and the diffu- 
sion terms are 



G [lm] ( P ) =n{ 7^(Qi)im h 



C [lm] ( P ) = 5R{ -w 



lm a lm 



}, 



(14) 
(15) 



D [ kKp) =®{ -uL • [/3D u (u)] im }, (16) 



D™(p) =5R{ - 7 C[^D e («)] im }. (17) 
The vertical fluxes and the surface term are 



-u)f m $i m - uf ro • [wu] im /2 }, (18) 



Fl™\p) =5R{ - 7 tM"V/2 }, 



(19) 



S 



+ul m ■ [® s V h p s 8] lm }, (20) 



where 6 = S(p s — p) is the Dirac distribution with impulse 
at the surface (see appendix). 

2.4 Vertical integration, summation over zonal wavenum- 
bers and cumulative summation over total wavenum- 
bers 

Since the density strongly varies with height in the atmo- 
sphere we prefer to include the density when we integrate 
spectral energies over layers, so that our quantities have 
the dimension of energy rather than energy per unit mass 
as in most other studies. With the formulation by Boer, 
this can be done easily even for pressure levels pierced by 
the topography. The vertically integrated KE spectrum is 
defined as 



EkW p \ 



rdp E 

J Vt 9 l<m<l 



(21) 
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The vertically integrated nonlinear spectral flux of kinetic 
energy is defined as 

n*=Ef - E r Jr ] (p). ( 22 ) 

and the spectral flux of APE is defined in the corresponding 
way. When © and pOJ) are vertically integrated, divided 
by g and summed as in (f2"2")l over all the spherical harmonics 
with total wavenumber ~> I, we obtain 



d t £ A m 



cm 



n f 



•^R:t['](P6 

cm 

FAf[l](Pb] 



em 



F Kt [l}( Pt )-V K {l}%{23) 
JU t [i](Pt)-^[e, (24) 



where the terms named J-"-j-[Z](p) = X)n>z -^tW are cumula- 
tive vertical fluxes and each of the other terms is an inte- 
grated cumulation of a corresponding term in Q and (fTUl) . 
For example, = Yl n >l is the cumulative 

kinetic energy and C/f[/]^ = J2 n>l C[n]p b t the cumulative 
conversion of APE into KE. Putting I = in equations 
(j!?3"j) and d25) we recover equations §5§ and (@| integrated 
between two pressure surfaces. 



3 Application to two GCMs 

3.1 Presentation of the data and the models 

We have analyzed two data sets produced by two sim- 
ulations performed with two spectral GCMs: the At- 
mospheric GCM for the Earth Simulator (AFES) and 
ECMWF's weather prediction model Integrated Forecast 
System ( IFS). For a prec i se de script ions of both sim ula- 
tions, see lHamilton et al.l (|2008h and IeCMWfI (|2010h . re- 
spectively. The two simulations and the two models are 
quite different. The AFES is a climate model using a 
spectral advection scheme. It has been run at high res- 
olution for research purposes. The horizontal resolution is 
T639, which correspond to a grid space of roughly 20 km. 
There are 24 vertical le vels from the ground to about 1 hPa. 
Hamilton et al. (|2008h showed that the AFES reproduces 



many features of the atmospheric spectra, especially a re- 
alistic fc -5 / 3 power law at the mesoscales. ECMWF IFS is 
a model developed and used for operational deterministic 
forecast. It uses a semi-Lagrangian advection scheme with 
a horizontal resolution T1279, which correspond to a grid 
space of roughly 10 km. There are 91 vertical levels with a 
minimum pressure of 0.01 Pa. The AFES and the ECMWF 
simulations correspond to June and December, respectively. 
We have averaged over 25 days (5 different times) for the 
ECMWF data set and over 10 days (40 times) for the AFES 
data set. 

The AFES data contain the horizontal velocity, the 
pressure velocity and the temperature at pressure levels. 
The geopotential is computed by integrating the hydro- 
static relation from the ground. The ECMWF data contain 
the vorticity, the divergence, the pressure velocity and the 
temperature at hybrid vertical levels. We compute horizon- 
tal velocity and the geopotential and then linearly interpo- 
late the data at pressure levels. Since we do not have access 



to the heating rate and to the total dissipative terms, the 
APE forcing G[l]p b t and the dissipative terms VK[l]p b t and 
T>A[l]p* have not been computed. Since the surface fluxes 
are modeled in a GCM, we focus on the vertical flux at 
pressure levels not pierced by the topography. Finally, we 
have computed the terms C[l]%, Ii- K [l] p p h t , U A [l]^ and C[l]f t 
for all levels and the vertical fluxes Fk\ [I] (p) and Fa\ [I] (p) 
for pressure levels not pierced by the topography. The spec- 
tral flux arising from the Coriolis and pressure terms £[l]p^ 
is completely negligible at wavenu mbers I > 30 consis- 
tent w ith previous computations bv lKoshvk and HamiltonI 
(2001). At larger scales, it is not negligible at one instant 
but shows a big randomness in time because it is a sum 
of large terms which tend to cancel each other. Therefore, 
one would need to average over long time series to obtain a 
good statistical convergence. Since our main interest here is 
the dynamics of the mesoscales and for clarity, this spectral 
flux is not included in the figures. 



3.2 Vertically integrated spectral energy budget 

Figure HJ^a) presents the globally integrated spectral fluxes 
and cumulative conversion for the AFES model. By con- 
struction the fluxes are equal to zero at I = Z max and should 
also be equal to zero at I — 0, since they represent conser- 
vative processes. The total spectral flux (black thick line) is 
positive at all wavenumbers, which means that in average, 
the energy is transferred toward large wavenumbers. At 
leading order, energy is forced at the very large planetary 
scales and dissipated at smaller scales, and a substantial 
part is dissipated at the smallest scales simulated. 

However, the total spectral flux has a somewhat intri- 
cate shape. It reaches a maximum equal to 1.3 W/m 2 
around I = 4, decreases to 0.55 W/m 2 at I ~ 20, increases 
again to reach a plateau between I ~ 70 and I ~ 200 where 
U[l] ~ 0.82 W/m 2 before dropping down to zero at the 
largest wavenumbers. At I < 15, the total flux is largely 
dominated by the APE spectral flux (blue line), which also 
increases abruptly around / = 3 and decreases abruptly be- 
tween I = 6 and I = 20. This indicates that there is a 
transfer of APE from wavelengths of the order of 10000 km 
to wavelengths between 2000 km and 5000 km. The strong 
decrease of [I] is associated with a strong increase of the 
KE spectral flux (red line) and a strong decrease of the cu- 
mulative conversion (dashed dotted line) from C[l] = 1.13 
W/m 2 to C[l] = -0.14 W/m 2 . The variation AC[AZ] of the 
cumulative conversion over a particular range of wavenum- 
bers represents the energy converted over this range. In 
this case, the conversion of APE to KE is AC[AZ] ~ 1.27 
W/m 2 over a range of wavelengths between 1000 km to 
5000 km. This large conversion at the synoptic scales and 
the spectral transfer of APE from the planetary scales to- 
ward the synoptic scales are mainly due to the baroclinic 
instability. It is interesting to compare these results with 
the spectral energy budget of an e quilibrated Eady flow 
( Molemaker and McWilliamsl . 2010h . The general picture 
is very similar with a dominance of the APE flux at large 
scales and a stronger KE flux at small scales. The increase 
of the total nonlinear flux at wavelengths between 700 km 
to 2000 km indicates that there is a direct forcing of APE 
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(a) wavelength (km) (b) wavelength (km) 




10° 10 1 10 2 10 3 '10° 10 1 10 2 10 3 

l + l l + l 



Figure 1: Total, KE and APE nonlinear spectral fluxes and cumulative conversion C[l] versus total wavenumber for (a) 
the AFES T639 simulation and (b) the ECMWF IFS T1279 simulation. The cross in (a) marks the maximum value of 
the KE nonlinear spectral flux Wk = max(IIx) = 0.62 used for nondimcnsionalization of the spectra in figure [21 In (b), 
the black dashed line is the spectral flux n[Z] for the AFES T639 simulation. 



at these scales, most probably due to latent heat release 
organized at la rge scales. This interp retation is consistent 
with results bv lHamilton et al . (2008). who reported spec- 
tral magnitude at the mesoscales much smaller for the dry 
dynamical core than for the full AFES. 

As already shown, the KE is mainly forced (by a con- 
version of APE) over a range of wavelengths between 1000 
km to 5000 km. At larger scales, the KE spectral flux 
is negative and reaches a minimum value approximately 
equal to -0.45 W/m 2 . A portion of the KE is transferred 
upscale, toward the planetary scales, which feeds the large- 
scale zonal winds. At smaller scales, the KE spectral flux 
reaches a plateau at 0.62 W/m 2 . This shows that a non- 
negligible portion of the KE cascades toward small scales 
at the mesoscales. The value 0.62 W/m 2 , which after con- 
version gives 6.1 x 10~ 5 m 2 /s 3 , is consistent with previous 
estimations for t he KE spectral flux and for the small-scale 
dissipation rate ( Cho and Lindborel 2001 ). Remarkably, 
the cumulative conversion C[l] increases at the mesoscales 
(the local conversion is negative), showing that the conver- 
sion is from KE to APE. This demonstrates that the KE 
£-5/3 gpectrum is not produced by direct forcing of the 
KE. Note that this conversion from KE to APE is consis- 
tent with strongly stratified turbulence. 

One should not interpret the wavenumber where IIr- [/] 
changes sign as the wavenumber where both upscale and 
downscale fluxes start. These two processes actually coex- 
ist in the same range of scales. The nonlinear KE spectral 
flux computed only with the rotational flow, Tl ro t, is plot- 
ted as a red dotted line. The interactions between the rota- 
tional modes conserve both KE and enstrophy (| V/j Au| 2 /2), 
exactly as in two-dimensional turbulence. Note that IL ro t 
is the spectral flux computed when the framework based 
on the two-dim ensional vorticity equation is adopted, as 



In contrast, the complementary flux, n^-[7] — Tl ro t, (red 
dashed line) is responsible for the downscale energy flux 
which starts at wavelengths of the order of 2000 km. In- 
deed, the downscale energy cascade is produced by inter- 
actions involving the d i vergent part of the v elocity field 
(jMolemaker et all l20iri iDeusebio et~aH . l2012h . 



As shown by iLambertl (|1984T ). the variations of the cu- 
mulative conversion at very small wavenumber I < 8 are 
due to the Hadley and Ferrel circulations. The decrease of 
C[l] at I = 2 corresponding to a conversion of APE to KE 
of approximately 0.6 W/m 2 is mainly the signature of the 
Hadley cell. The increase of C[l] at I = 4 corresponding 
to a conversion of KE to APE of approximately 0.5 W/m 2 
is mainly the signature of the Ferrel cell. The c onversion 
at I = 2 is weaker than previously computed by lLambert 
. Further investigations are necessary to understand 



for example in iBoer and Shepherd! (|1983l ). We see that 
these interactions are responsible for the upscale flux, which 
actually starts at wavelengths of the order of 1800 km. 



if this effect is due to averaging over an insufficient amount 
of data or if it is a robust aspect of the AFES. 

Figure [TJb) presents the globally integrated spectral 
fluxes and cumulative conversion as in figure Ufa) but for 
the ECMWF model. The large-scale features are overall 
quite similar to the results for the other model even though 
the conversion at wavenumbers smaller than 4 correspond- 
ing to the Hadley cell is much stronger (AC[A^] ~ 1 W/m 2 ). 
The baroclinic instability leads to large-scale positive APE 
flux and strong conversion at the synoptic scales AC[AZ] ~ 
1.5 W/m 2 . However, the KE flux increases much less than 
for the AFES, indicating that there is strong dissipation 
at wavelengths of the order of 1000 km. The total spec- 
tral flux for the AFES is also plotted in dashed black line 
for comparison. The total flux for the ECMWF model is 
slightly smaller than for the AFES at the synoptic scales 
and is negative and very small at the mesoscales. This is 
related to the weakness of the downscale mesoscale KE cas- 
cade (U K [l] ^0.15 W/m 2 ) and to the fact that the APE 
flux is negative (n^Z] ~ —0.3 W/m 2 ). This unexpected 
result could be due to direct forcing of APE by release of 
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latent heat at the smallest resolved scales of the order of 
50 km. This interpretation is consistent with the sign of 
the local conversion at the mesoscales, from APE to KE, 
in contrast to the case of the AFES. 

3.3 Vertically integrated non-dimensional spectra 

If it is assumed that the Z~ 5 / 3 range of the kinetic energy 
spectrum is of a similar form as the Kolmogorov spectrum 
of isotropic turbulence, then Ek oc n 2 ^ 3 / -5 / 3 . If it is fur- 
ther assumed that the only other parameters that deter- 
mine the spectrum are the Earth's radius, a, and the total 
mass per unit area, which is equal to (p s )/g, then dimen- 
sional considerations give 



E K [i) = c({ Ps )/g) 1 ^(an K ) 2/3 i- 5/3 , 



(25) 



where C is a constant supposedly of the order of unity. In 
the following, we choose as the typical KE flux the maxi- 
mum of the KE flux in the mesoscales: Hk = max(IlA') = 
0.62. (This value is marked by a cross in figure 
In figure dJ^a) the non-dimensional compensated spectra 
E[l}l b ' 3 ((p s ) / g)~ 1 / 3 (atl K )- 2/:i for the AFES model are 
plotted as a function of the total wavenumber. At / = 1, the 
APE spectrum (blue line) is muc h larger than t he KE spec- 
trum (red line), as predicted bv iLorenz ( 19551 ). This leads 
to a ratio mean APE over mean KE approximately equal to 
3. However, for other wavenumbers (except at the largest 
ones), both spectra are of the same order of magnitude. 
At the synoptic scales, the KE spectrum is quite steep. It 
shallows at the mesoscales and presents a flat plateau cor- 
responding to a l~ 5 / 3 inertial range from 650 km up to 
the large-wavenumber dissipative range. Remarkably, the 
constant C in equation (|25J) is very close to unity, which 
indicates that the /~ 5 / 3 mesoscale ran ge may be explained 
in a similar way as iKolmogorovl ( 1941 ) explained the fc~ 5 / 3 
range of isotropic turbulence. At wavelengths between 700 
km to 2000 km, the APE spectrum (blue line) is equal to 
or larger than the KE spectrum and there are fluctuations 
resembling noise. However, the fluctuations do not seem re- 
lated to lack of statistics. They could be due to the direct 
forcing of APE at this scales (see figure [T^.). 

In figure [5] are also plotted the rotational spectrum 
E rot [l] and the divergent spectrum Edi v [l]- In order to de- 
fine these spectra, one has to consider an alternative ex- 
pansion of the velocity field than in ([5]), 

U = E E ^^(xd+^Yg^W (26) 



l>0 -l<m<l 



where Y 



(rot) 



ae z A V h Y lm and Y% v) 



aV h Yi m are the 



vector spherical harmonics satisfying the orthogonal con- 
dition (Y^yJ,^,) = 1(1 + l)5 aa ,5w5 mm >. The spherical 
stream function tj}(yih,p) is related to the vertical compo- 
nent of the vorticity by £ = e z ■ (V/i Au) = —a\Wh\ 2 4' an d 
the spherical velocity potential xi^hiP) is related to the 
horizontal divergence by d = Vft • u = a|V^| 2 x- We see 
that the divergent and rotational spectra can be defined as 



Ediv [I] 



E 



a 2 \d h 



21(1+1)' 



E rn f \l\ 



^ 21(1 + 1)' 



With these definitions, J2i>o Ediv[l] + E rot [l] is equal to 
the total energy although EdivW + E rot \l] is not exactl y 
equal to the KE spectrum ( Koshvk and Hamilton! . l200lh . 
At large scales, the rotational spectrum E rot [l] (dashed 
red line) largely dominates over the divergent spectrum 
Ediv[l] (dashed dotted red line). Remarkably, the com- 
pensated divergence spectrum increases with Z, meaning 
that Ediv[l] is shallower than a l~ 5 / 3 power law. Such 
very shallow divergent spectra were also obtained by sim- 
ulations of strongly stratified and st rongly rotating turbu- 
lence forced in geostrophic modes ( Deusebio et all [2012) 



and of s trongly stratified tu rbulence forced with rotational 



modes ( Augier et al. . 2012f h However, other simulations 
with higher resolution would be necessary in order to 
check whether the divergent spectra are sensitive to model 
parametrizations and resolution. At the mesoscales, both 
spectra are of the same order of magnitude even though 
E rot [I] is still larger than Edi V [I] ■ 

Figure [2Jb) shows the non-dimensional compensated 
spectra for the ECMWF model (we use the same value 
as for the AFES of fl^ in order to allow a easier compar- 
ison between both models). The large-scale features are 
very similar to the AFES. At the synoptic scales, the ra- 
tio Ek [1}/Ea[1] is approximately equal to 2, indicating that 
the energy is partitioned equally between the two compo - 
nents of KE and the APE, as predicted bv lCharnevI (1971). 
In contrast to the KE spectrum, the APE spectrum be- 
comes more shallow at high wavenumbers. This is prob- 
ably related to direct forcing of APE at the mesoscales. 
Interestingly, the compensated divergent spectrum is flat, 
which means that it follows a l~ 5 / 3 power law. However, 
its magnitude is very small so that the vertically integrated 
KE spectrum at the mesoscales is nearly not affected. It is 
interesting to note t he similarity with t he fc~ 5 / 3 divergent 
spectra obtained bv lWaite and Snvderl ( 20091 ) simulating a 
baroclinic life cycle. 



3.4 Vertical decomposition and vertical energy fluxes 

Figure [3] presents the spectral fluxes and the cumulative 
conversion integrated over two different layers correspond- 
ing approximately to the upper troposphere (a,b) and the 
stratosphere (c,d). Figures (a,c) correspond to the AFES 
T639 simulation and figures (b,d) to the ECMWF IFS 
T1279 simulation. Figure [3J also presents the cumulative 
total vertical fluxes, J-f[l](pb), at the bottom, and J-^[Z](pt), 
at the top of the layer (magenta dashed and dashed-dotted 
lines, respectively). The balance between these two inward 
and outward terms, A^'J^fZ] = J-f[l] (pb) — ^[l] (pt), is plot- 
ted in magenta as a continuous line. 

For the AFES simulation, the spectral fluxes and the cu- 
mulative conversion integrated over the upper troposphere 
(figure present roughly the same features as the glob- 
ally integrated terms shown in figure [IJa). The KE flux 
in the upper troposphere accounts for approximately one 
third of the globally integrated KE spectral flux. The cu- 
mulative vertical flux A^J 7 ^] is relatively small at the 
mesoscales and at the synoptic scales, indicating that at 
leading order, the energy budget in this layer and at these 
scales is dominated by the spectral fluxes rather than by 
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Figure 2: Non-dimensional compensated spectra versus total wavenumber for (a) the AFES T639 simulation and (b) the 
ECMWF IFS T1279 simulation. The black dashed line represents the prediction (j2l>|) with Tlx = 0.62 (value marked by 
a cross in figure [TJl) . The continuous straight line indicates the l~ 3 power law. 




Figure 3: Same as figures [T] but integrated over layers corresponding approximately to (a,b) the stratosphere and (c,d) the 
upper troposphere. Figures (a,c) correspond to the AFES T639 simulation and figures (b,d) to the ECMWF IFS T1279 
simulation. The legend is given in (a) and A^J^[Z] = J-f[l](j>b) — J~\[l]{pt) is the cumulative inward vertical flux. 
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the vertical fluxes. However, at the top and the bottom 
of the layer, there are large vertical fluxes which are ap- 
proximately equal. These fluxes through the upper tropo- 
sphere (upward at wavenumbers 7 $5 I ^ 20 and downward 
at wavenumbers 3 ^ I ^ 6) account for exchanges of en- 
ergy between the lower troposphere and the stratosphere. 
At wavenumbers 1 ^ I ^ 2, F^[l](pt) decreases from 1 
W/m 2 to W/m 2 indicating a strong upward vertical flux 
at pt = 233 hPa. Since the vertical flux at the bottom layer 
is small at these wavenumbers, the layer loses energy. The 
wavenumber ranges of the upward and downward fluxes at 
the planetary scales nearly coincide with the wavenumber 
ranges of the conversion due to the Hadley and Ferrel cells, 
which seems to indicate that these vertical fluxes are related 
to the large-scale cells. In contrast, the upward flux at the 
synoptic scales indicates the presence of upward propagat- 
ing planetary waves. At I > 70, F^[l]{jpb) and T^[l](pt) 
increase meaning that there are downward vertical fluxes. 
This shows that in the AFES, the mesoscales of the upper 
troposphere are not directly forced by upward propagat- 
ing gravity waves, as was t he case in the simulation by 



Koshvk and Hamilton! (|200lh using the SKYHI model 



Figure [3fb) presents the same quantities as figure O^a) , 
also integrated over the upper troposphere but for the 
ECMWF model. Comparing figures |3Ja) and[3Jb), we see 
the same differences as we saw between figures []Ja) an d 
[ljb). The vertical fluxes are also quite different from the 
AFES with small upward fluxes at the mesoscales and a 
smaller magnitude of the variations of the cumulative fluxes 
at large scales. 

Figure [2£c) shows the same quantities as figure G^a) , 
also for the AFES but integrated over the stratosphere. At 
the large scales, the spectral fluxes and the cumulative con- 
version are quite different from the terms integrated over 
the height of the atmosphere and over the upper tropo- 
sphere. The APE flux has no large peak at I ~ 4 and there 
is a conversion from KE to APE at synoptic scales rather 
than the other way around, as in the upper troposphere. 
On the other hand, there is a strong upward energy flux at 
large scales from the bottom layer at pb — 233 hPa while 
the corresponding flux at the top layer at p t = 15 hPa is 
very small. Thus, the stratosphere is not directly forced by 
baroclinic instability but by an energy flux from the tropo- 
sphere, due to upwards propagating planetary waves and 
possibly also the effects of the Hadley and the Ferrel cells 
(see figure [5k). At the mesoscales, there is a conversion 
of KE into APE and a strong downscale cascade of KE 
and APE accounting for approximately half the globally 
integrated flux. At I > 80, A^J^[Z] ~ ^[^(pb) increases, 
meaning that the stratosphere loses energy by a downward 
flux through the tropopause. 

Figure ED(d) presents the same quantities as figure [3Jc), 
also integrated over the stratosphere but for the ECMWF 
model. Note that the vertical axis is different from fig- 
ure E3c). In contrast to the AFES, there is nearly no down- 
scale energy cascade and the stratospheric mesoscales are 
directly forced by an upward energy flux. 
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E K , p > 412 hPa 
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E K , p < 233 hPa 

E A , p < 233 hPa 




Figure 4: KE and APE spectra integrated over three lay- 
ers corresponding approximately to the lower troposphere, 
the upper troposphere and the upper atmosphere. The 
straight black lines indicate the and l~ 5 ^ 3 dependencies 
as guides for the eye. Note that the axes and the straight 
lines are exactly the same in (a) and (b). 



Figure U] presents the KE and APE spectra integrated 
over three layers corresponding approximately to the lower 
troposphere, the upper troposphere and the upper atmo- 
sphere, for the AFES (figure |Sl) and the ECMWF model 
(figure 0Jd) . Note that the axes and the straight lines are 
exactly the same in both figures. At the large scales the 
spectra for both models are similar but at the mesoscales 
the spectral magnitude of the ECMWF is much smaller 
than the magnitude of the AFES. Spectra integrated over 
the different layers lie closer to each other for the AFES 
than for the ECMWF. For the AFES the KE spectra have 
similar shapes at different layers, whereas for the ECMWF 
model, they are quite different with a shallowing at the 
mesoscales only in t he stratospher e (blu e and red dotted 
lines) as reported bv lBurgess et al. ( 20121) . However, these 
shallow spectra are much smaller in magnitude than the 
corresponding spectra from the AFES model (more than 
one order of magnitude smaller at I ~ 100). From fig- 
ure [BJd) it is quite clear that the shallowing of the strato- 
spheric spectra of the ECMWF model is not caused by 
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an energy cascade but rather by gravity waves propagat- 
ing from the troposphere. The integrated spectra of the 
ECMWF model (figure [5b) is dominated by the contribu- 
tion from the troposphere. In the lower troposphere, the 
APE spectrum presents a l~ 5 ^ 3 dependency and is larger 
than the KE spectrum at large wavenumbers. However, this 
cannot be explained by a downscale energy cascade since 
the APE spectral flux is small and negative. These differ- 
ences between the APE and the KE spectra could be due 
to a damping of the KE and to a direct diabatic forcing of 
the APE at the very large wavenumbers I ~ 1000. For the 
AFES model (figure 2k), the spiky irregular shape observed 
at the synoptic scales for the globally integrated APE spec- 
trum (figure Hk) is also observed for the tropospheric APE 
spectra but not for the stratospheric APE spectrum. This 
could indicate that this effect is related to the topography. 
However, the spikes are not only present in the lower but 
also in the upper tropospheric spectra whereas the upper 
troposphere is not pierced by the topography. Moreover, 
the vertical fluxes at p — 412 hPa are rather regular. There- 
fore, the spiky irregular shape of the APE spectra seems to 
be due to processes happening in the troposphere and not 
only at the surface. The release of latent heat organized at 
the synoptic scales could be an explanation. 

4 Conclusions 



velocity field (IVallgren et all 1201 It iDeusebio et all 12012 : 
Molemaker et all l201oI ~ one must consider the nonlinear 
transfers computed with the total advection term. To do 
this in a consistent way, the pressure term and the advec- 
tion terms in the spectral energy budget have to be decom- 
posed as shown in section [5J 

Future investigations could compare different models 
varying the resolution, the convection schemes, the advec- 
tion scheme and the turbulent models. It would be de- 
sirable to analyze simulations over long periods to obtain 
better statistical convergence at large scales. In particular, 
this would be necessary in order obtain statistical conver- 
gence of the linear term ()13|> at large scales. Moreover, it 
could be informative to analyze idealized simulations us- 
ing dry dynamical core and/or aquaplanet versions of the 
models with variation of physical parameters such as the 
rotation rate and the large-scale forcing. Other interesting 
aspects for future work are the consideration of the effects 
of the water content and the adaptation of the formulation 
for non- hydrostatic simulations. 

We wish to thank Kevin Hamilton and Nils Wedi for 
providing the data. We also thank Kevin Hamilton for 
showing nice hospitality in Hawaii and Nathanael Schaeffer 
for pro viding his l ibrary for spherical harmonic transform 
SHTns <|Schaeffeill20l5i . 



A new formulation of the spectral energy budget has been 
presented and applied to study the results of two differ- 
ent GCMs. In contrast to previous formulations, both KE 
and APE are considered and the topography is taken into 
account. Moreover, the advection terms are exactly sepa- 
rated into spectral transfer and vertical flux and the pres- 
sure term is exactly separated into adiabatic conversion, 
vertical flux and spectral transfer. 

The spectral fluxes show that the AFES, which re- 
produces realistic spectra with a k~ 5 ^ 3 power law at the 
mesoscales, simulates a strong downscale energy cascade of 
II[Z] ~ 0.8 W/m 2 . The vertical fluxes for the upper tropo- 
sphere show that the mesoscales are not directly energized 
by gravity waves propagating from the ground. Moreover, 
the spectra collapse on the prediction based on the exis- 
tence of the cascade, indicating that the mesoscale fc~ 5 / 3 
power law is due to the downscale energy cascade. In 
contrast, neither the fc _5//3 spectra nor the downscale en- 
ergy cascade are produced by the ECMWF model. The 
study of the spectra and their tendencies integrated over 
different layers reveals that in the AFES, the stratospheric 
mesoscales are not forced by upward propagating gravity 
waves. In contrast to this, the stratospheric spectra of the 
ECMWF model are directly forced at the mesoscales by 
gravity waves propagating from the troposphere. 

These results show that our spectral energy budget for- 
mulation is a convenient tool to investigate the issue of the 
mesoscale dynamics and its simulation by GCMs. In partic- 
ular, we have shown that the flux computed only with the 
rotat i onal part of the velocity fie ld n rot (as i.e. in lFiOrtoftl . 
19531 : Boer and Shepherd . 1983) accounts only for the up- 
scale KE flux. Since the downscale energy cascade is pro- 
duced by interactions involving the divergent part of the 



Technical details on topography in 
p-coordinates 



Following iBoerl (1982), the equations are extended over a 
domain including subterranean pressure levels. For any 
function f(xh,p) whose values below the surface have 
been obtained by interpolation, we define a correspond- 
ing function f(x h ,p) = /3(xfo,p)/(x h ,p), where P(xh,P) = 
H(p s (xfj) — p) and H is the Heaviside function. In the ex- 
tended domain, the boundary condition at the earth surface 
is expressed as D t( 3 = 0, with D f = dt + v • V. Using the 
chain rules d t (3 = Sd t p s , = SVhPs, d p (3 — —S, we re- 
cover the classical boundary condition for the atmospheric 
domain in pressure coordinates: D t f3 — 5(dtp s + u • Vh.p s — 
lo s ) = 0, where 8 = 5{p s —p) is the Dirac distribution with 
impulse at the surface. 

In practice, computing the spherical harmonics trans- 
form of u = /3u is not numerically feasible at pressure lev- 
els pierced by the topography, i.e. where j3 is equal both 
to 1 and to 0. We must use a modified smooth /3. More- 
over, the spherical harmonic transform of u = /3u should 
not reflect the spectral content of the topography rather 
than that of u. Therefore, it is convenient to use an al- 
ternative /3 computed w ith a time-averaged pressure sur- 



face p s (yih) (|Boerl . 119821 ). A smooth low-pass filter with 



a cutoff total wavenumber equal to 40 is then applied to 
this function H(p7(x/ l ) — p). In this way, a large part of 
the interpolated subterranean data (mainly under the large 
topographic highs, i.e. Antarctic continent, Himalaya and 
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Andean mountain ranges) are not used in the calculations, 
but the spectral quantities at relatively high wavenumbers 
are not affected by the high wavenumber content of the 
topography. 
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